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A  continuum-level  electrochemical  model  previously  developed  by  the  authors  [1  ]  is  used  to  investigate 
the  dependence  of  open-circuit  voltage  (OCV),  and  maximum  power  density  on  electrolyte  thickness  for 
solid  oxide  fuel  cells  (SOFCs)  with  mixed  conducting  electrolytes.  Experimental  results  confirm  the  mod¬ 
els  predictions  that  OCV  decreases  monotonically  with  decreasing  electrolyte  thickness  due  to  increased 
permeation  flux  [1].  The  model  was  further  extended  to  show  that  there  exists  an  optimal  electrolyte 
thickness  at  which  maximum  power  density  occurs  for  mixed  conducting  electrolytes.  As  expected,  for 
electrolyte  thickness  greater  than  optimal  losses  from  ohmic  overpotential  reduce  cell  output.  How¬ 
ever,  when  the  electrolyte  thickness  is  lower  than  optimal  losses  from  an  increasing  electronic  “leakage” 
current  reduce  cell  output. 
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1.  Introduction 

Acceptor-doped  ceria  has  widely  been  investigated  as  an  elec¬ 
trolyte  in  solid  oxide  fuel  cells  (SOFCs),  especially  at  intermediate 
temperatures  (500-700  °C)  [1-7].  However,  although  acceptor- 
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doped  ceria  has  excellent  ionic  conductivity  (up  to  5  times 
higher  than  yttria-stabilized  zirconia,  YSZ  [5-7])  it  is  a  mixed 
conductor— i.e.,  it  conducts  electrons  in  addition  to  oxygen  ions— in 
typical  SOFC  operating  conditions  [5,6].  This  results  in  a  delete¬ 
rious  leakage  current  during  operation  that  reduces  overall  SOFC 
efficiency  [2-6].  This  leakage  current,  as  a  fraction  of  total  current, 
increases  as  the  cell  voltage  increases  and  in  open-circuit  condi¬ 
tions  its  effect  shows  up  as  a  reduction  in  the  open-circuit  voltage 
from  the  theoretical  (Nernst)  value. 
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Recent  studies  have  shown  that  electrolyte  thickness  is  also 
a  factor  in  the  magnitude  of  the  leakage  current  and,  therefore, 
the  open  circuit  voltage  [1,8].  Indeed,  it  was  shown  experimen¬ 
tally  [8]  and  theoretically  [1]  that  reducing  electrolyte  thickness 
also  reduced  open-circuit  voltage.  However,  open-circuit  voltage 
is  a  marker  for  the  maximum  potential  energy  available  to  an 
SOFC.  Hence  for  SOFCs  with  mixed  conducting  electrolytes  the 
typical  approach  of  making  ever  thinner  electrolytes  to  increase 
maximum  current  density  (by  reducing  ohmic  losses)  may  not  be 
applicable  to  MIEC  electrolytes.  In  other  words,  for  mixed  conduct¬ 
ing  electrolytes  there  is  likely  a  trade-off  between  reducing  ohmic 
losses  and  maintaining  a  high  enough  open-circuit  potential.  Con¬ 
sequently,  one  may  reasonably  infer  that  there  exists  a  optimum 
thickness  at  which  maximum  power  density  is  attained.  SOFC  elec¬ 
trolytes  with  a  thickness  below  the  optimum  value  have  a  reduced 
power  density  because  of  too  great  a  loss  in  available  potential 
energy,  i.e.,  too  low  an  open-circuit  voltage.  Similarly,  SOFC  elec¬ 
trolytes  with  a  thickness  above  the  optimum  value  also  have  a 
reduced  power  density,  but  now  because  of  increased  ohmic  losses. 

In  addition,  a  survey  of  the  literature  shows  great  inconsistency 
in  the  open-circuit  voltages  (OCVs)  reported  by  various  groups 
for  SOFCs  with  acceptor-doped  ceria  electrolytes,  even  with  elec¬ 
trolytes  of  similar  thickness  [1-7].  While  the  inconsistency  in  OCVs 
could  be  attributable  in  part  to  differences  in  GDC  film  quality  (due 
to  variations  in  porosity,  homogeneity,  etc.),  they  are  nevertheless, 
particularly  stark  when  contrasted  to  the  absence  of  similar  vari¬ 
ation  for  SOFCs  with  YSZ  electrolytes  (which  have  similar  issues 
with  film  quality).  Since  YSZ  is  a  purely  ionic  conductor,  one  may 
deduce  that  the  variability  seen  in  SOFCs  with  acceptor-doped  ceria 
electrolytes  is  also  a  consequence  of  mixed  conductivity  therein. 

As  a  corollary  of  the  authors’  previous  work  on  modeling  of  SOFC 
electrochemical  performance  [1],  this  paper  explores  the  hypoth¬ 
esis  introduced  above  by  modeling  the  performance  of  real  SOFCs 
from  our  laboratory  and  elsewhere  in  an  effort  to  deduce  the  opti¬ 
mal  thickness  of  acceptor-doped  ceria  electrolytes  and  to  properly 
describe  all  the  factors  that  play  a  role  in  measured  open-circuit 
voltages  of  SOFCs  with  acceptor-doped  ceria  electrolytes. 

2.  Theory 

A  general  model  for  SOFC  performance  in  all  operating  condi¬ 
tions  was  derived  in  entirety  in  our  earlier  work  [1  ].  Modeling  the 
dependence  of  open-circuit  voltage,  @0c ,  on  thickness  and  electrode 
characteristics  is  a  subset  of  that  work  and  so,  in  this  section,  only 
the  relevant  features  are  summarized. 

2.1.  Current-voltage  relationships 

The  equivalent  circuit  of  an  SOFC,  of  thickness  L,  is  shown  in 
Fig.  1,  where  oxygen  partial  pressure,  Po2,  o  (at  x  =  0,  i.e.,  the  anode 
side)  and  Po2,l  (at  x  =  L,  i.e.,  the  cathode  side)  are  the  Po2 ’s  on  either 
side  of  the  MIEC.  Fig.  1  shows  the  most  widely  reported  SOFC  con¬ 
ventional  equivalent  circuit  in  the  literature  [9-16].  Circuit  analysis 
shows,  given  p  =  po  +  t)l,  that 

i  =  Jv  +Je  =  -  r,  -  <Pext Wu  '  -  (1 ) 

where  p  is  the  total  interfacial  activation  overpotential,  pQ  is  the 
activation  overpotential  at  the  anode-electrolyte  interface,  pL  is  the 
activation  overpotential  at  the  cathode-electrolyte  interface,  J  is 
current  density  drawn  by  the  load  (i.e.,  external  circuit),  @ext  is  the 
potential  difference  across  a  load  or  applied  by  an  external  source, 
0th  is  the  theoretical/Nernst  voltage  and  R  is  electrical  resistance. 
Additionally,  the  subscripts  “0”  and  “L”  refer  to  the  anode/reducing 
side  (i.e.,  at  x  =  0)  and  the  cathode/oxidizing  side  (i.e.,  at  x  =  L)  of  the 
electrolyte,  respectively  (as  shown  in  Fig.  1);  while  the  subscripts 
“ v ”  and  “e”  refer  to  oxygen  vacancies  (i.e.,  the  mobile  ionic  species) 


x  =  0 


x  =  L 


Fig.  1 .  Equivalent  circuit  of  an  SOFC.  0th  represents  the  theoretical  or  Nernst  voltage. 
The  cell  output  voltage,  is  applied  to  a  load  resistance,  RL,  through  which  flows 
the  cell  output  current  J,  electronic  ‘leakage’  current  Je,  flows  through  the  electronic 
resistance,  Re,  while  the  cell  activation  overpotentional  is  represented  by  Ract,o  and 
Ract,L  for  the  anode  and  cathode,  respectively.  The  ionic  current  Jv  is  greater  than  J 
and  flows  through  the  ionic  resistance  Rv. 


and  electrons,  respectively.  Fig.  1  and  Eq.  (1 )  both  show  that  at  open 
circuit  C/  =  0)JV  =Je  ±  0  and  &oc  =  (®th  -rj-  (Pext)ReR^  Therefore, 
since  Re  and  Rv  are  both  functions  of  thickness,  it  is  clear  that  @oc 
will  be  as  well. 

2.2.  Thermodynamics 

For  dilute  concentrations  of  point  defects  (i.e.,  where  the  elec¬ 
trochemical  potential  of  lattice  oxygen  jl0x  is  approximately 

constant),  the  oxygen  exchange  reaction  between  the  MIEC  and 
the  gas  phase  is  given  (in  Kroger-Vink  notation)  by 

Oq  ~  V"  +  2ft  +  \o2,  K  =  cvc2epfj  (2) 

where  K  is  the  mass  action  (equilibrium)  constant  for  the  reaction 
and  c  is  concentration. 

At  the  anode  interface,  the  reaction  in  Eq.  (2)  proceeds  to  the 
right,  while  at  the  cathode  interface  it  proceeds  to  the  left.  Accord¬ 
ingly,  it  was  shown  in  our  earlier  work  [1  ]  that 

I<0  =  1C  exp(-2qkg1T_1  no)  =  P) 

I(L  =  K*  exp(2qkg1T_1  m)  =  c^cfpff  (4) 

where  I<0  and  I<L  are  the  equilibrium  constant  for  oxygen  exchange 
at  the  anode-electrolyte  interface  and  cathode-electrolyte  inter¬ 
faces,  respectively,  and  is  the  equilibrium  constant  for  I<*  oxygen 
exchange  when  po  =  pi  =  0\  also  T  is  temperature,  q  is  elementary 
electron  charge  and  1<b  is  Boltzmann’s  constant.  Eqs.  (3)  and  (4) 
show  that  I<  is  affected  by  the  overpotential  at  the  electrode,  i.e.,  the 
amount  of  oxygen  incorporated  (cathode)  in  or  released  (anode)  by 
the  electrolyte  is  a  function  of  the  electrode  overpotential.  Conse¬ 
quently,  through  Eq.  (2),  the  concentration  of  defect  species  at  the 
boundaries  (i.e.,  at  x  =  0  and  x  =  L)  is  dependent  on  potential.  It  is 
noted  here  that  defect  formation  in  electrode  materials  is  excluded 
because  it  primarily  affects  electrode  kinetics  which  is  discussed  in 
the  following  section. 

2.3.  Kinetics 

The  transport  model  derived  in  the  next  section  must  be  related 
to  the  kinetics  of  electron  transfer  at  the  electrodes.  This  is  done 
through  the  Butler-Volmer  equation,  which  may  be  approximated 
as  follows  [17,18] 

jv  =  2Jex  0  sinh(2_1q/<:“1T_177o)  relative  to  the  anode  (5) 

or 

jv  =  2/ex, i  sinh(2_1q/<“1T_1r/L)  relative  to  the  cathode  (6) 
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where  Jex  is  the  exchange  current  density.  In  the  Butler-Volmer 
equation,  Jex  encapsulates  the  material  properties  of  an  electrode 
that  determine  its  performance  in  the  SOFC. 

Moreover,  if  a  is  defined  as  the  ratio  of  the  anodic  exchange 
current  density  to  the  cathodic  exchange  current  density  (i.e., 
a=Jex,olJex,L)>  then  combining  Eqs.  (5)  and  (6)  yields 


2.5.  Open  circuit  voltage 

In  open  circuit  conditions,  there  is  no  current  drawn  by  an  exter¬ 
nal  circuit  so  that  J  =  0.  From  the  authors’  earlier  work  [1],  the 
following  expressions  for  <P0C  and  the  open-circuit  oxygen  vacancy 
current,  JVt0C,  were  derived, 


r]L  =  2q  U<bT  sinh  1  (a  sinh  Qqkg’r  %))  (7) 

i.e.,  the  cathodic  overpotential  is  a  function  of  the  anodic  overpo¬ 
tential  and  vice  versa. 

This  approach  is  a  significant  departure  from  existing  models 
[9-16]  in  which  the  distinctions  between  the  individual  charac¬ 
teristics  (as  determined  by  their  material  properties)  of  the  anode 
versus  the  cathode  are  ignored.  Specifically,  the  Butler-Volmer 
equation  is  employed  in  such  a  way  that  only  p  is  reported  (i.e., 
the  distinct  contributions  of  p0  and  pL  are  ignored)  and  Jex  is 
representative  of  the  effective  exchange  current  density  of  both 
electrodes  combined.  Arguably,  this  practice  may  be  acceptable  for 
ideally  reversible  electrodes,  but  for  real,  irreversible  electrodes 
more  sophistication  is  desirable.  In  the  proposed  approach  out¬ 
lined  above,  it  is  possible,  in  principle,  to  separate  p  into  po  and 
Pi ,  and  obtain  Jex0  and  JeXyL.  We  hasten  to  point  out,  however,  that 
this  approach  can  only  provide  apparent  values  for  p0  and  pi ,  and 
thereby  JeXy0  and  JeXyL ,  since  the  absolute  measurement  of  p0  and  pL 
requires  the  use  of  a  suitable  reference  electrode  [17-20].  Never¬ 
theless,  the  apparent  values  obtained  for  p0  and  pL ,  and  thereby ]eXi o 
and  JeXtL,  will  provide  much  needed  insight  into  the  performance 
characteristics  of  SOFCs. 


®oc  =  &th  - 


I<bT  f  1  cVjo  De  —  Dv  zv(ue  —  Uv)Cv,o  —  ueca  \  M 
1100  q  U  cv,L  +  zvDv  -  zeDe  n  zv(ue  -  uv)cv,l  -  ueca  )  ' 


L  = 


uvueqXy0Cca 
2{uv-ue]lex,0  . 
qXyocuvueCa 


sinh(2 b^00) 


Jv’°C~  ( uv-ue)L 
given 

A  =  fcsTq-V 


=  2/ex,o  sinh  (^^o.oc^ 


and 


(11) 

(12) 


Yoc  —  Y  ( ^oc  +  Poc  —  & th  —  ~Z  ln  )  —  (cv,L  —  Cy, o)  (13) 
A  V  zvq  cVy0 


where  L  is  electrolyte  thickness,  z  is  charge  number  (ze  =  -l  and 
z„  =  2),  D  is  diffusivity,  u  is  mobility  and  ca  is  acceptor  dopant 
concentration— preset  by  material  synthesis  and  for  the  case  study 
material,  ceria  (Ce02 )  the  acceptor  dopant  of  choice  is  trivalent  Gd3+ 
or  Sm3+. 

Eqs.  (10)  and  ( 1 1 )  are  a  set  of  parameterized  equations  in  terms 
of  p0yOC.  Flence  the  dependence  of  @oc  on  L  can  be  educed  by  first 
obtaining  the  O0c{po,oc)  and  L(po,oc)-  This  was  done  with  JeXy0  and 
JeXyL  as  the  only  fitting  parameters. 


2.4.  Physical  significance 

The  crucial  feature  of  the  thermodynamic  and  kinetic  equations 
developed  above  is  the  occurrence  of  p  in  each  of  the  equation 
sets:  Eqs.  (3)  and  (4)  (thermodynamics)  and  Eqs.  (5)-(7)  (kinetics). 
For  mixed-conducting  electrolytes  this  is  fundamental,  because  it 
shows  how  electrode  properties  influence  the  thermodynamics 
(i.e.,  the  defect  populations)  of  the  electrolyte  (and  vice  versa).  This 
becomes  clear  after  substituting  Eqs.  (5)  and  (6)  into  Eqs.  (3)  and 
(4),  respectively,  as  shown  below 

-1  (\  Jv_ 

y  2Jex,0 

I<L  =  1C  exp  (4  sinh-1  A. J ^  =  cv,LclLPf22L  (9) 

Consequently,  as  will  demonstrated  later,  this  approach  allows 
the  characteristics  of  any  type  of  electrode  to  be  correctly  incorpo¬ 
rated  into  the  SOFC  performance  model  developed  herein. 

The  effect  is  clearly  seen  in  <P0c  measurements  reported  for 
SOFCs  with  mixed  conducting  electrolytes  such  as  acceptor-doped 
ceria.  In  open-circuit  conditions  with  mixed  conducting  elec¬ 
trolytes,^  0  but  Jv  =f=  0  hence  p  =f=  0.  Therefore,  as  Eqs.  (8)  and  (9) 
show,  the  defect  populations  (as  well  as  I<0  and  I<L )  are  dependent 
on  Jex> o  and  JeXiL  even  at  open-circuit.  In  other  words,  with  a  mixed 
conducting  electrolyte,  @0c  depends  on  the  electrodes  employed. 
This  result  is  supported  by  the  literature  [21-23]  where  (for  the 
same  T,  Pq2,o  and  Pq2  L )  there  is  an  appreciable  variance  in  0OC 
reported  for  SOFCs  with  acceptor-doped  ceria  electrolytes  which 
is  largely  due  to  the  electrodes  used,  Fig.  1.  In  contrast,  for  SOFCs 
with  YSZ  electrolytes  the  variance  in  reported  @oc  is  negligible 
since,  for  purely  ionic  conducting  electrolytes  like  YSZ ,J=JV  =  0  and 

p  =  0. 


2.6.  Maximum  power  density 

Electrolyte  thickness  also  plays  a  role  in  maximum  power  den¬ 
sity.  In  pure  ionic  conducting  electrolyte,  the  maximum  power 
density  increases  monotonically  with  decreasing  electrolyte  thick¬ 
ness.  Flowever,  for  mixed  conducting  electrolytes  (due  to  the  higher 
mobility  of  the  electronic  species),  it  is  possible  that  the  magni¬ 
tude  of  the  leakage  current  increases  with  decreasing  thickness 
faster  than  the  ionic  current.  Therefore,  with  mixed  conducting 
electrolytes,  the  SOFC  output  power  can  end  up  eventually  decreas¬ 
ing  as  the  electrolyte  thickness  decreases.  This  phenomenon  can  be 
explored  by  plotting  the  maximum  power  density  as  a  function  of 
thickness  using  the  equation  below  [1]: 

Pload  —  ^extJ  =J  ( th  ~  Pact  ~  Relec  +  ~~z  J 

\  ZVQ  Cy, 0  J 

|  zv(De  -  Dy)y  +  (J/q)  ln  Zv{ue  -  Uy)cVyL  -  ueca  -  (JL/qXy ) 
zv(uv  -  Ue){y/JL)  zv{ue  -  uv)cVy0  -  ueca  -  {JL/qXy ) 

(14) 

where  Reiec  is  the  sum  of  the  non-electrolyte  ohmic  resistances, 
including  the  resistances  of  the  leads,  the  current  collectors  and 
the  electrodes  as  well  as  contact  resistances;  &*xt  is  therefore  the 
corrected  external  potential. 

As  detailed  by  Duncan  and  Wachsman  [  1  ] ,  Eq.  ( 1 4)  can  be  used  to 
obtain  the  current-voltage  and  current-power  curves  a  fuel  cell  of  a 
given  electrolyte  thickness.  Consequently,  after  fitting  the  model  to 
experimental  data  for  a  given  electrolyte  thickness,  the  maximum 
power  density  can  be  predicted  for  other  electrolyte  thicknesses 
(assuming  concentration  overpotential  is  negligible). 

3.  Experimental 

To  validate  the  models  described  above  a  series  of  anode  sup¬ 
ported  SOFCs  were  fabricated  with  varying  GDC  thickness.  Special 


K0  =  IC  exp  (  -4  sinh 


=  cv,ocl  0P{ 


o2,o 


(8) 
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Fig.  2.  SEMs  of  the  cell  cross-sections  for  electrolyte  thicknesses  of  (a)  5.3  |jim,  (b)  9.4  |jim,  (c)  11.3  |jim,  (d)  12.6  (Jim  and  (e)  32.9  (Jim.  In  each  figure,  the  top  layer  is  LSCF-GDC 
cathode,  the  bottom  layer  is  the  Ni-GDC  anode  and  the  middle  layer  is  the  electrolyte.  The  anode  functional  layer  (AFL)  is  not  visible  at  these  magnifications. 


care  was  taken  to  fabricate  the  cells  with  high  repeatability  and 
steps  were  taken  to  minimize  the  variability  between  cells.  Specifi¬ 
cally,  the  anode  supports  were  obtained  from  the  same  tape,  while 
the  anode  functional  layer  (AFL)  and  cathode  inks  used  were  taken 
from  the  same  batch. 

3.1.  Cell  fabrication 

NiO-GDC  anode  supports  were  prepared  by  tape-casting;  using 
a  mixture  of  NiO  (Alfa  Aesar)  and  GDC  (Rhodia  powder)  with  a  65 :35 
weight  ratio.  NiO  and  GDC  powders  were  ball  milled  with  Solsperse 
as  a  dispersant  in  a  mixed  toluene/ethyl  alcohol  solvent  system  for 
24  h.  Next,  a  plasticizer  and  a  binder  were  added  to  the  suspension 
and  ball  milled  for  another  24  h.  A  mixture  of  di-n-butyl  phthalate 
(DBP)  and  polyethylene  glycol  (PEG),  and  polyvinyl  butyral  (PVB) 
were  used  as  plasticizer,  and  binder,  respectively.  The  slurry  was 
then  degassed  in  a  vacuum  chamber  to  prevent  the  formation  of 
cracks  or  defects  from  air  bubbles  during  the  tape-casting  process. 
Then  a  continuous  NiO-GDC  anode  support  tape  was  fabricated 
from  the  slurry  on  a  Procast®  tape  casting  system  (DHI,  Inc).  After 
drying  for  2  h  at  1 00  °C,  the  NiO-GDC  tape  was  punched  into  circular 
green  bodies  with  32  mm  diameter  and  partially  sintered  at  900  °C 
for  2  h  to  get  rid  of  residual  organic  binders  and  provide  proper 
mechanical  strength  for  the  electrolyte  deposition  process. 

NiO-GDC  colloidal  slurry,  containing  submicron  size  NiO  (JT 
Baker®)  and  GDC  (Rhodia  powder)  (65:35  wt  %),  was  mixed  with 
organic  binders  and  coated  on  one  side  of  anode  support  by  spin 
coating  to  form  an  anode  functional  layer  (AFL)  for  the  SOFC.  (The 


spin  coater  used  was  set  at  1500rpm  for  15  s.)  Subsequently,  the 
anode  substrates  were  pre-heated  at  400  °C  for  the  removal  of  the 
binders. 

To  fabricate  dense  GDC  electrolytes,  a  GDC  colloidal  slurry  was 
prepared  and  deposited  by  spin  coating  Rhodia®  GDC  powder  was 
ball  milled  for  24  h  with  Solsperse®  (dispersant)  in  ethanol.  PVB 
(binder)  and  DBP  (plasticizer)  were  added  after  the  first  ball-milling 
step  and  the  solution  was  ball-milled  for  an  additional  24  h.  The  spin 
coating  was  the  same  as  for  the  AFL  deposition.  In  order  to  con¬ 
trol  the  thickness  of  GDC  electrolyte,  the  spin  coating  process  was 
repeated  for  each  sample  until  the  desired  thickness  was  obtained. 
After  deposition,  samples  were  dried  at  room  temperature  for  10  h 
and  the  multilayered  fuel  cell  structure  was  sintered  at  1450  °C  for 
4  h  in  air. 

Lao.6Sro.4Coo.2Feo.803_5  (LSCF)-GDC  composite  cathode  were 
prepared  and  applied  on  the  GDC  electrolyte  surface.  Cathode  inks 
were  synthesized  by  mixing  LSCF  (Praxair® )  and  the  GDC  (Rhodia®) 
in  a  50:50  weight%  with  organic  binders.  After  mixing  and  grinding 
the  cathode  ink  for  1  h,  the  ink  was  brush-painted  onto  the  GDC 
electrolyte  evenly.  After  applying  the  ink  3  times,  the  cathode  was 
sintered  at  1100° C  for  1  h.  The  active  cathode  area  was  ~0.4cm2. 
For  current  collection,  Ag  mesh  and  Pt  wire  were  attached  onto  both 
electrode  surfaces  using  Pt  paste  and  then  fired  at  900  °C  for  1  h. 

3.2.  Electrochemical  characterization 

Samples  were  loaded  in  a  fuel  cell  testing  set-up  using  517- 
ceramabond  sealant  (Aremco®),  30  seem  of  dry  air  and  30  seem  of 
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Fig.  3.  Fit  of  the  OCV  model  (Equation)  to  experimental  data  for  open-circuit  voltage 
(OCV)  as  a  function  of  electrolyte  thickness  (I)  for  a  Ni-GDC/GDC/LSCF-GDC  cell. 


wet  hydrogen  were  supplied  to  the  cathode  and  anode  side,  respec¬ 
tively.  After  reduction  of  anode  for  several  hours,  the  open  circuit 
voltage  (OCV)  and  current-voltage  (/-V)  characteristics  at  various 
temperatures  were  obtained  using  a  Solartron®  1287  potentiostat. 

4.  Results  and  discussion 
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Fig.  4.  Arrhenius  plot  of  fitting  parameters  used  to  model  the  Ni-GDC/GDC/LSCF- 
GDC  cell. 


the  plots  show  that  if  the  electrolyte  is  made  sufficiently  thick  the- 
oretical/Nernst  value  for  OCV  is  achievable;  albeit  at  great  cost  to 
performance  in  closed-circuit  (]  ±  0)  conditions. 

The  continuum-level  electrochemical  model  for  SOFC  per¬ 
formance,  Eq.  (14),  was  fitted  to  experimental  data  for  a 
Ni-GDC/GDC/LSCF-GDC  SOFC  (with  an  anode  functional  layer)  in 
Fig.  6a.  Normally,  the  model  uses  three  fitting  parameters,  I<,  a  (i.e., 


Fig.  2  shows  the  SEMs  of  the  cell  cross-sections  for  various  elec¬ 
trolyte  thicknesses,  from  5.3  |xm  to  32.9  [xm.  In  each  figure,  the  top 
layer  is  LSCF-GDC  cathode,  the  bottom  layer  is  the  Ni-GDC  anode 
and  the  middle  layer  is  the  electrolyte.  The  anode  functional  layer 
(AFL)  is  not  visible  at  these  magnifications. 

The  OCV  model,  Eqs.  (10)  and  (11),  was  fitted  to  experimental 
data  of  open-circuit  voltage  as  a  function  of  electrolyte  thickness 
for  a  Ni-GDC/GDC/LSCF-GDC  cell,  for  three  temperatures  (500  °C, 
600  °C  and  650  °C),  Fig.  3.  The  plot  shows  the  model  approximates 
the  data  quite  well  given  the  error  inherent  in  the  experiment.  The 
results  clearly  show  that  as  the  electrolyte  thickness  decreases, 
so  does  the  OCV  for  the  cells.  Moreover,  the  decrease  in  OCV 
per  decrease  in  electrolyte  thickness  increases  as  L  decreases, 
i.e.,  AOCV/AL  ocl/L.  Physically,  this  result  suggests  an  exponen¬ 
tial  increase  in  the  permeation  flux  as  the  electrolyte  gets  thinner. 
And,  therefore,  the  reduction  in  ohmic  resistance  afforded  by  a  thin¬ 
ner  electrolyte  is  countered  by  a  decrease  in  convertible  energy  (as 
indicated  by  the  lower  OCV). 

The  two  parameters  (Jex<0  and  Jex<L )  extracted  from  fitting  the 
model,  were  plotted,  Fig.  4,  to  verify  Arrhenius  behavior,  since 
kinetic  phenomena  are  normally  parameterized  by  an  activation 
energy.  The  excellent  fits  obtained  from  the  plots  (correlation 
coefficient,  ft  >0.98)  confirm  Arrhenius  behavior,  but,  more  sig¬ 
nificantly,  they  confirm  the  consistency  of  the  OCV  model.  In 
practical  terms,  requiring  that  the  extracted  parameters  conform 
to  Arrhenius  behavior  increases  their  meaningfulness.  The  activa¬ 
tion  energies  for  the  cathode  and  anode  exchange  current  densities 
are  comparable  (1.1  eV  for  the  anode  and  0.94  eV  for  the  cathode) 
and  reasonable  for  these  materials  [24]. 

Fig.  5a  shows  the  oxygen  permeation  flux  through  the  GDC  elec¬ 
trolyte  (in  open-circuit,  i.e.,  zero  current  conditions)  and  OCV  as  a 
function  of  thickness,  at  650  °C.  Correspondingly,  Fig.  5b  shows  a 
plot  of  OCV  versus  the  permeation  flux,  at  650  °C.  From  the  plots  it 
is  easy  to  see  that  as  the  electrolyte  thickness  decreases  the  oxy¬ 
gen  permeation  flux  increases  while  the  OCV  decreases.  Moreover, 
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Fig.  5.  (a)  OCV  and  oxygen  permeation  flux  as  a  function  of  electrolyte  thickness  at 
650  °C;  (b)  OCV  as  a  function  of  oxygen  permeation  flux  at  650  °C. 
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Jex,olJex,L)>  and  Reiec  [1  ].  However,  a  K  value  of  4.64  x  1066  m9  atm1/2 
was  independently  obtained  from  thermogravimetry  experiments 
[25]  on  the  same  GDC  powder  used  in  the  electrolyte  and  the 
electrodes.  So  with  only  two  fitting  parameters  ( a  and  Reiec)  good 
agreement  between  model  and  data  was  obtained.  This  is  a  signifi¬ 
cant  achievement  considering  contemporary  performance  models 
all  use  more  (>5)  fitting  parameters  and  in  some  cases  fit  only  about 
half  the  performance  range  (i.e.,  from  open-circuit  to  maximum 
power  density)  [13,15,16]. 

The  model  returns  a  value  of  a  =  1 .1 9  which  means  the  effective 
exchange  current  densities  for  the  anode  is  slightly  better  than  the 
cathode  for  this  cell.  In  addition,  this  value  of  a  is  almost  identical 
to  the  value  used  to  fit  the  OCV  vs.  L  data,  at  600  °C,  in  Fig.  3,  which 
was  oi  =Jex,o IJex,L  =  1 .27.  The  difference  between  the  two  values  is 
6.3%,  which  is  well  within  the  boundaries  of  reasonable  experimen¬ 
tal  error.  In  addition,  the  model  yields  Reiec  =  0.001  £2  cm2,  which 
is  quite  a  small  value,  that  demonstrates  the  effect  of  the  AFL  on 
lowering  the  effective  ohmic  resistance  of  the  anode. 

The  parameters  used  in  fitting  the  Ni-GDC/GDC/LSCF-GDC  SOFC 
performance,  in  Fig.  6a,  were  then  used  to  predict  the  dependence 
of  the  maximum  power  density  (MPD)  of  that  cell  as  a  function  of 
electrolyte  thickness  in  Fig.  6b.  The  plot  shows  that  (at  MPD)  though 
the  ohmic  overpotential  (which  is  directly  dependent  on  electrolyte 
thickness)  decreases  monotonically  with  decreasing  electrolyte 
thickness,  the  electronic  leakage  current  increases  (exponentially) 
with  decreasing  electrolyte  thickness.  Consequently,  as  the  elec- 
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Fig.  6.  (a)  Model,  Eq.  (14),  fitted  to  experimental  data  for  current-voltage  and 
current-power  of  a  Ni-GDC/GDC/LSCF-GDC  SOFC  (with  an  anode  functional  layer) 
with  a  12-p.m  thick  electrolyte  at  600  °C.  (b)  Projected  maximum  power  density, 
electronic  (leakage)  current  and  ohmic  overpotential  as  a  function  of  electrolyte 
thickness  predicted  from  model  at  600  °C. 
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Fig.  7.  Projected  maximum  power  density,  electronic  (leakage)  current  and  ohmic 
overpotential  as  a  function  of  electrolyte  thickness  predicted  from  model  at  600  °C 
for  (a)  Ni-SDC/SDC/BSCF  [21]  and  (b)  Ni-GDC/GDC/LSCF-GDC  [22]  cells. 


trolyte  thickness  decreases  the  MPD  initially  increases  but  when 
the  leakage  current  becomes  too  great  the  MPD  reverses  course  and 
begins  decreasing,  after  passing  through  a  local  maximum.  In  other 
words,  there  exists  an  optimal  electrolyte  thickness  for  a  mixed 
conducting  electrolyte  such  as  GDC.  This  concept  is  crucial,  because 
it  means  that  mixed  conducting  electrolytes  should  not  be  made 
arbitrarily  thin.  Reducing  the  electrolyte  thickness  increases  the 
electronic  leakage  current  and  decreases  the  OCV,  both  of  which 
mitigate  and  can  even  overwhelm  the  benefit  of  reduced  ohmic 
overpotential. 

Fig.  7  shows  the  MPD  dependence  on  thickness  generated  for 
previously  modeled  SOFCs  [21,22].  The  basic  trends  obtained  in 
Fig.  6b  are  also  seen  for  these  SOFCs.  However,  the  shape  of  the 
curves  were  influenced  by  the  electrode  materials  used  for  those 
SOFCs,  which  therefore  have  different  values  for  a  Jo, ex  >  and  I<.  These 
differences  were  discussed  thoroughly  in  our  earlier  work  [1  ]. 


5.  Conclusions 


The  continuum-level  electrochemical  model  for  SOFC  perfor¬ 
mance  was  applied  to  the  dependence  of  open-circuit  potential  and 
maximum  power  density  on  electrolyte  thickness  for  SOFCs  with 
acceptor-doped  ceria  electrolytes. 

Good  fits  of  the  model  were  obtained  for  experimental  data  for 
OCV  versus  thickness  for  different  temperatures.  Furthermore,  the 
results  showed  that  OCV  decreases  with  decreasing  L.  This  means 
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that  for  mixed  conducting  electrolytes  the  ideal  SOFC  electrolyte 
thickness  is  a  compromise  between  the  need  to  lower  ohmic  losses 
(thinner  L)  and  the  desire  to  maximize  convertible  energy  (higher 
OCV).  The  fitting  parameters  were  meaningful  and  conformed  to 
expected  Arrhenius  behavior. 

The  continuum-level  electrochemical  model  was  fitted  to  SOFC 
performance  (current-voltage,  current-power)  data  using  only 
two  fitting  parameters  (down  from  three)  because  of  the  availabil¬ 
ity  of  accurate  independent  data  for  I<  from  thermogravimetry.  This 
is  a  significant  improvement  over  competing  models  that  use  more 
fitting  parameters  ( >5)  thereby  reducing  the  physical  and  statistical 
significance  of  each. 

The  SOFC  performance  parameters  were  then  used  in  the  model 
to  predict  the  dependence  of  the  maximum  power  density  as  a 
function  of  thickness  for  various  SOFCs.  The  findings  reveal  that 
there  is  an  optimal  thickness  at  which  the  power  density  is  the 
greatest.  For  thicknesses  less  than  optimal  the  power  density 
decreases  due  to  high  leakage  current  densities.  Conversely,  for 
thicknesses  higher  than  optimal  the  power  density  is  reduced  due 
to  growing  ohmic  losses. 

In  summary,  the  presence  of  electronic  conductivity  in  mixed 
conducting  electrolytes  means  that  their  thicknesses  should  not  be 
chosen  to  be  arbitrarily  thin.  Instead,  their  thickness  must  be  part 
of  the  design  strategy  to  maximize  the  SOFC  performance  require¬ 
ments. 
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